A quantum computation architecture using optical tweezers 
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We present a complete architecture for scalable quantum computation with ultracold atoms in 
optical lattices using optical tweezers focused to the size of a lattice spacing. We discuss three 
different two-qubit gates based on local collisional interactions. The gates between arbitrary qubits 
require the transport of atoms to neighboring sites. We numerically optimize the non-adiabatic 
transport of the atoms through the lattice and the intensity ramps of the optical tweezer in order 
to maximize the gate fidelities. We find overall gate times of a few 100 /is, while keeping the error 
probability due to vibrational excitations and spontaneous scattering below 10 -3 . The requirements 
on the positioning error and intensity noise of the optical tweezer and the magnetic field stability are 
analyzed and we show that atoms in optical lattices could meet the requirements for fault-tolerant 
scalable quantum computing. 



I. Introduction 

Neutral atoms in periodic optical potentials have long 
been considered a promising candidate for scalable quan- 
tum computation due to long coherence times of inter- 
nal state qubits and their excellent controllability. How- 
ever, a great challenge in these systems remains to tailor 
interactions for the two-qubit gates. One approach is 
to use dipole-dipole interactions 1J, most prominently 
among Rydberg atoms [2], realized recently with indi- 
vidual atoms in separate dipole traps Hj. A differ- 
ent approach uses ground state collisions either within 
the same trapping potential 5 - 8 or mediated by tunnel- 
ing [5] [TU]. In optical lattices, two-qubit gates have been 
implemented on many pairs of atoms in parallel |lll 112] 
but they have, so far, not been realized on a single pair 
of atoms. The main challenge is to attain sufficient reso- 
lution to manipulate single lattice sites. Several propos- 
als to overcome the diffraction limit have been made to 
achieve this goal [T5lll6| . and recently single-site manip- 
ulation by an optical tweezer in a short-period optical 
lattice has been realized [T7] and could now be exploited 
for the implementation of single- and two-qubit gates. 

In this work, we propose and investigate an architec- 
ture for scalable quantum computation in optical lat- 
tices. Using tightly focussed optical tweezers, atomic 
qubits are transported around in the lattice and merged 
into single lattice sites in order to implement collisional 
quantum gates. We adapt the ^/swap gate, proposed 
in Ref. [18] and realized in Ref. [12] to our tweezer-based 
optical-lattice architecture. By moving the atoms non- 
adiabatically, our transport and two-bit gate times be- 
come much faster than, e.g., the tunneling gate proposed 
in Ref. [S]. Similar to work in Ref. [T5], we numerically 
optimize the ramp-up of the tweezer intensity and the 
transport, but due to the difference in the available con- 
trol parameters, we obtain smaller errors, and arrive at 
total gate times of a few hundred microseconds with er- 



ror probability 10~ 3 arising from non-adiabatic excita- 
tions and spontaneous scattering. This scheme can in 
principle allow for thousands of operations within the 
coherence time of the qubit. We also investigate the re- 
quirements for the stability of the position and intensity 
of the tweezer as well as the effect of magnetic field noise. 

The paper is organized as follows. Section II intro- 
duces the scalable architecture, including the quantum 
register initialization, the single-qubit gates and the read- 
out. Section III investigates the non-adiabatic transport 
of atoms through the lattice, which is required for the 
collisional two-qubit gates. In Section IV, we present 
three different gates and discuss their respective advan- 
tages including the role of coherence times of the qubit 
states used. In section V, the different errors introduced 
by spontaneous scattering and by fluctuations of tweezer 
position and intensity are discussed and compared to the 
threshold for fault-tolerant quantum computing. 
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Figure 1: Optical tweezer-based architecture, a) A tightly 
focused optical tweezer locally perturbs the optical lattice po- 
tential, b) The atoms can be shuttled around by moving the 
tweezer, c) The quantum register is initialized with an atomic 
pattern with empty rows and columns to permit unhindered 
transport of atoms within the lattice. The two-qubit gates 
utilize the collisional interactions between atoms positioned 
at the same site. With multiple tweezers in use, several gates 
can be performed in parallel. 



II. Scalable architecture 

In this section, we present our scalable architecture for 
quantum computing using optical tweezers for address- 
ing and transport in optical lattices. The differential light 
shift of a tweezer on selected atoms can be exploited for 
single-qubit gates (Section II. C), while moving tweezers 
can bring arbitrary pairs of atoms together for collisional 
two-qubit gates. Single-qubit rotations form a univer- 
sal set of gates in combination with either the two-qubit 
phase gate pj] or the ^/swap gate [20]. Together with 
the initialization (Section II. B) and the readout capabil- 
ity (Section II. D), the system therefore fulfills the require- 
ments for scalable quantum computation 21j. 



A. General setup 

Throughout this paper we assume a 3D optical lat- 
tice of 87 Rb atoms formed by three pairs of counter- 
propagating near-infrared laser beams (wavelength Ai a t — 
1064 nm), realizing a lattice with spacing ai a t = Ai a t/2. 
We assume the ability to prepare a single 2D system and 
the availability of an imaging system with a resolution 
such that all individual sites can be detected [22] [23] 
and individually addressed [T7] [53]. We will assume a 
lattice depth in all directions of V\ a t = 50 E T , where 
E r = h 2 /(2mAiat) is the recoil energy and m is the 
mass of a Rb atom. This yields a trapping frequency 
w trap ~ 27T • 30 kHz, a two-particle interaction energy 
Ui n t/h s=s 2 kHz, and a tunnel coupling J/h ps 0.06 Hz. 

The qubit will be composed of two ground state hy- 
pcrfine levels of 87 Rb but unlike in Ref. [TTJ we will as- 
sume the atoms to be addressed by a tweezer close to 
the 5S-6P transition at X 1 / 2 = 421.555 nm and A 3 / 2 = 



420.1733 nm. Tuning the wavelength between these two 
transitions yields a differential light shift which we will 
exploit for spin-dependent transport. We describe the 
optical tweezer by a Gaussian beam profile with a 1/e 2 
waist radius of w = ai a t/2, which can be reached with 
an imaging system of numerical aperture 0.68. Unless 
otherwise stated the maximum light shift of the optical 
tweezer is V t — 500 E r . For the gate in Section IV. C, two 
tweezers are needed, and many tweezers are required for 
the parallelization of gates. These can be created and 
quickly controlled using acousto-optic deflectors [25] . 



B. Quantum register initialization 

A major advantage of ultracold atoms in optical lat- 
tices over other systems lies in the possibility of creat- 
ing large scale quantum registers. This can be done by 
preparing the system in a Mott insulating state (551 [23J 
[25] , which pins the occupation of the lattice sites to inte- 
ger values due to the on-site interaction between atoms 



and can thus realize unity occupation. When the lattice 
is sufficiently deep, deviations from unity filling are only 
thermally activated and depend on the position in the 
external confinement. In the center of the Mott insulat- 
ing domain, the deviations can be very small and current 
experiments reach P n j=\ ~ 3T0~ 3 [25] ■ Different purifica- 
tion schemes have been proposed to circumvent residual 
thermal defects. They involve either illumination and 
subsequent recooling of all atoms [27J , many-body dy- 
namics [2BTl31| . or algorithmic cooling schemes [32]. Re- 
cent experiments have demonstrated a number-selective 
removal of atoms with an efficiency of 86% [33J . 

Starting from the Mott insulator with unity filling, we 
remove every second column and additionally one row so 
that they can be used as a channel for shuttling atoms 
between different sites [Fig. [It)]. The removal can be 
performed by transferring selected atoms to a different 
hyperfine ground state and subsequently removing them 
with a resonant laser illuminating the entire sample |17) . 



C. Single-qubit gates 

A single-qubit gate can be realized using either mi- 
crowave radiation or optical two-photon Raman transi- 
tions between two selected ground state hyperfine levels. 
In the former case, the differential light shift of a tightly 
focused tweezer or a magnetic field gradient shifts se- 
lected atoms into resonance with a global microwave field 
driving the transition. In the latter case, the rotation is 
driven by two tightly focussed laser beams in a Raman 
configuration with zero two-photon detuning. 

Raman-based gates have been demonstrated in sin- 
gle optical dipole traps and reached a 7r/2 gate time of 
183 ns [33] and 37 ns [35] with a next-neighbor residual 
rotation at the level of 10 -3 [34]. A theoretical analy- 
sis for the performance in a large spacing optical lattice 
yields an achievable error of 10~ 4 for Rb [35] and 10 -5 
for Cs [37]. 

Microwave-based gates were realized in a ID lattice us- 
ing a magnetic field gradient, reaching a 7r/2 gate time 
of 8/xs 14 1 . In 2D, single spin manipulation was demon- 
strated using the local light shift due to a focused laser 
beam and Landau-Zener microwave sweeps [17] , A theo- 
retical analysis suggests that gate times of 100 /is can be 
reached with an error probability of ending in the wrong 
state below 10~ 4 at a lattice spacing of 10 /im [37] and 
with an error < 10~ 3 at a spacing of 425nm[24. Re- 
cently an error of 1.4 • 10~ 4 per global microwave-based 
single-qubit gate of 30 /xs gate time was demonstrated 
for Rb atoms in an optical lattice [35] using randomized 
benchmarking |39j . 

Comparing the approaches, it is obvious that the 
Raman-based gates are much faster. However, since even 
the slower microwave-based gates are comparable with 
the two-qubit gate time discussed in Section IV faster 



single-qubit rotations will not significantly speed up the 
overall computation. The drawback of the Raman-based 
gate is that is extremely difficult to avoid cross-talk at 
the small lattice spacings envisaged here. On the other 
hand, the microwave-based gate restricts the choice of the 
qubit states to magnetic field-sensitive states, which have 
shorter coherence times, but which are anyway required 
for some two-qubit gates, such as the one discussed in 
Section IV.B. 



D. Readout 

As the final step, the spin state of several qubits has to 
be read out. For this, we can transfer one of the two qubit 
spin states into the neighboring free shuttling areas using 
the spin-dependent transport discussed in Section IV.B, 
followed by spatially resolved fluorescence imaging [23l 
HUI WT\ . In the case of magnetic field- insensitive qubit 
states, one can remove one of the spin states and detect 
the qubit state by the presence or absence of an atom in 
the fluorescence image. 



III. Non-adiabatic transport of atoms 

Our proposal is based on the collisional interaction of 
atoms brought to the same lattice site, and high-fidelity 
transport of atoms through the lattice is therefore cru- 
cial. The transport of a single atom to an empty lattice 
site is illustrated in Fig.!!]}). The optical tweezer is first 
focussed on the site containing the atom and its power is 
ramped up. Then its position is translated to the destina- 
tion site and the power is ramped down. Since collisional 
gates are sensitive to the vibrational state, it is crucial 
to avoid vibrational excitation in each step. In the fol- 
lowing we first review the free-space harmonic oscillator 
theory and subsequently numerically model the process 
in the presence of the lattice potential. We show that 
dynamics during the ramp-up is essentially unchanged, 
whereas the fidelity of non-adiabatic transport is signif- 
icantly reduced by the presence of the lattice potential. 
We then numerically optimize the parameters controlling 
the transport and thereby regain high fidelity transport. 



A. Optical tweezer ramp-up 

In the harmonic oscillator approximation, the vibra- 
tional excitations during ramp-up of the optical tweezer 
power can be treated analytically |52]. Because the ma- 
trix elements of the corresponding Hamiltonian vanish 
between states with opposite parity, the change of the 
trapping frequency can only induce excitations to the 
second excited state. In the low-excitation limit, the 
time-dependent Schrodinger equation is then solved in 
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Figure 2: Excitation probability versus the total ramp-up 
time T r of an optical tweezer centered on a lattice site to 
a depth of 500 E T according to Eq. (flj> . The inset shows the 
trapping frequency of the potential well as a function of time. 



the Hilbert space consisting of only the ground state \4> g ) 
and the second excited state \4> e ), separated by the en- 
ergy gap AE ge = 2Tluj. 

The adiabaticity criterion for a change in the trap fre- 
quency is 



dcu(t) 



dt 
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(AE g 



dH 
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where 



dH 



\4> g )\ = h/^/2 with the Hamiltonian H of 
an harmonic oscillator with varying trap frequency ui{i). 
If we keep the adiabaticity factor £ -c 1 constant, the 
above differential equation determines the shape of the 
power ramp-up. The total ramp-up time T T from an ini- 
tial trap frequency uj to a final trap frequency ujf is 
related to the adibaticity factor via 021 
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and the excitation probability is given by 
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The excitation probability displays an oscillatory behav- 
ior with an envelope P c ° = 4£ 2 /(l +4£ 2 ). The oscillatory 
factor in Eq. (13J) predicts that with an appropriate tim- 
ing, the ramp-up can be done non-adiabatically with in 
principle unity fidelity. The existence of minima in the 
non-adiabatic excitation probabilities have been demon- 
strated experimentally using a cloud of atoms in a dipole 
trap with a longitudinally translated focus [33"] . 

Since realistic lattice potential wells are not fully har- 
monic, we have studied the ramp-up of a Gaussian- 



shaped optical tweezer potential superimposed on a co- 
sine lattice potential. We numerically solve the time- 
dependent Schrodinger equation as the depth of the 
tweezer is increased from zero to 500 E T according to 
the adiabaticity condition of Eq. ([!]). 

The resulting excitation probability Pg Ull (T r ) as a func- 
tion of the total ramp-up time T T (solid line) decreases 
for larger ramp-up times illustrating the transition from 
non-adiabatic to adiabatic dynamics (Fig. [2]). For com- 
parison, we also show the envelope of the harmonic oscil- 
lator approximation P ° with £ calculated from the given 
total ramp time T r via Eq. (pi) (dashed line in Fig. [2]) and 
observe a decent but not exact agreement with the nu- 
merical result. The small disagreement is due to the an- 
harmonicity of the lattice potential and we have verified 
that it diminishes as the optical tweezer is made tighter. 
Our calculation also shows that very small excitation er- 
rors can be reached with carefully chosen ramp timing. 
For example, errors below P,!" 11 = 10~ 3 are realized at 
the first minimum at a ramp- up time of T x — 11 /is. 



B. Atom transport 

Atoms are transported in the lattice by shifting the po- 
sition x of the optical tweezer. Again, we are interested 
in the excitation probability during this process and we 
first discuss the transport of atoms in a harmonic po- 
tential without the presence of the lattice, which can be 
solved analytically. The effects of non-adiabaticity can 
be calculated by replacing w with x in Eq. (fTl) and using 
shifts the parity of states. In the 



^Q•- 



AE ge — hui, as 



8H 

c)x,, 



harmonic approximation 



; |H|0 9 > = hu/(y/2v ), 

where u = y/h/(mui) is the harmonic oscillator length. 
Again, we start with a translation profile which keeps 
the adiabaticity parameter £ constant. This requires 
a displacement with a constant velocity v 
£v2cr a-> and yields the excitation probability 



dt 



P e harm (i) = P c ° • sin 2 y/\ + 4£ 2 w£/2 



(4) 



with P ° as defined above. To estimate the time for an 
adiabatic transport over one lattice site, we consider the 
adiabaticity parameter £ = 0.016 for which the envelope 
P° of the excitation probability has dropped to 10~ 3 . 
For a potential depth of 500 E T of the optical tweezer 
(^trap ~ 2-7T • 90 kHz, <j — 36 nm), this yields a velocity 
of v — 0.45 /im/ms or a transport time of T t — vq ■ ai a t — 
1.2 ms over the distance of a single lattice spacing. 

We now include the lattice and the Gaussian profile of 
the optical tweezer and investigate to which extent non- 
adiabatic transport can be realized. Fig. pa) shows our 
numerical solution of the excitation probability Pg UU for 
a translation of the tweezer position over one lattice site 
at a constant speed as a function of the total transport 
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Figure 3: Excitation probability versus transport time, a) 
Linear displacement of the optical tweezer. Numerical result 
-Pe" 11 for the excitation probability including the influence of 
the lattice and the profile of the tweezer (solid line) and enve- 
lope P,? of the result for the transport of an harmonic oscil- 
lator potential without a lattice (dashed line), b) Zoom into 
the second minimum in a) . Numerical result for the linear dis- 
placement (dashed line) and an optimized ramp shape (solid 
line). The optimization reduces the excitation probability to 
below 10 -4 . The inset shows the displacement curve of the 
linear ramp (dotted line) and the optimized ramp (solid line) 
including harmonics up to the fifth order. 



time T t . Again, we observe an oscillatory behavior and a 
good agreement of the upper envelope curves. However, 
the numerical results do not reach zero excitation prob- 
ability during the oscillations. An enlarged view of one 
of the minima at short transport times shows that the 
excitation probability stays quite high (Fig.^)). 

To improve this, we parameterize the translation pro- 
file as a linear ramp plus a number of harmonics and 
optimize their weights using a standard simplex algo- 
rithm [44 j . Using only five harmonics dramatically im- 
proves the fidelity [solid curve in Fig.^)]. As can be 
seen in the inset, the translation profile is only slightly 
modified compared to the linear profile and the dominant 
component is a single harmonic with a period correspond- 



ing to the duration of the displacement. 

We have investigated how the error accumulates as the 
qubit is transported for longer distances across multiple 
sites. We have simply repeated the optimal single-site 
displacement found above with T t — 25 /is per lattice 
site up to 100 times. We find that the error after multi- 
site transport does not accumulate with the number of 
sites. Instead, it oscillates with a period of approximately 
eight lattice sites and is bound by a maximal error of 
only a few times the error of a transport over a single 
lattice site. Thus, fast shutteling of atoms within the 
quantum register should be feasible with total excitation 
probabilities well below 10~ 3 . 

For the transport tweezer we would choose a wave- 
length red detuned compared to both excited state levels 
(Dl and D2 line) such that both qubit states are trans- 
ported simultaneously If the polarization is chosen to be 
linear, both states experience the same light shift and do 
not undergo a differential phase evolution (the residual 
differential light shift drops below 10~ 3 for a detuning 
of A = 4 nm) . The effect of the transport beam on the 
neighboring sites can be neglected [24] . 



IV. Two-qubit gates 

In the section above we have described how we can po- 
sition two selected atoms in neighboring sites as a start- 
ing point for a two-qubit gate. The gate can be accom- 
plished by collisional interactions between atoms in the 
same well, and we will analyze three distinct proposals 
based on controlled atomic collisions, utilizing either the 
spin-dependent interaction, spin-dependent transport or 
spin-exchange interactions. 

The quantum bit is encoded in a pair of atomic states, 
one from each hypcrfinc split ground state manifold in 
87 Rb. We denote the upper and lower states |f) and 
|4), and in the following two subsections we show how a 
pair of atoms can acquire a phase evolution depending 
on the states of the two qubits, \jk) — > e %u ik l l h \jk), for 
j, k =t, I, and how the interaction strengths Ujk can be 
suitably controlled. 



A. Merging two qubits into a single well 

For the first two-qubit gate we discuss, one merges the 
two atoms into the ground state of a single combined 
well. The atoms have the same wave function and the 
interaction energy is calculated from an integral over the 
fourth power of the wave function. For cold atoms, the 
short range interaction strength is proportional to the s- 
wave scattering length, and near a Feshbach resonance 
there can be a difference between the scattering lengths 
for different and identical spin states. For the states 
\F =l,m F = +l) = |t) and |F = 2,m F = -l) = \i) 



of Rb, the interaction strengths at a magnetic field of 
B = 9.12 G are given by U n « U u « 0.9 U n 45J. Thus, 
by timing the merging and separation of the wells one can 
realize a phase gate. For an optical tweezer of potential 
depth 500-E r , as considered here, a differential phase of 7r 
between the two hyperfine states is acquired in roughly 
0.8 ms. 

A major experimental limitation with this approach is 
that during the merging and splitting of the atoms, the 
atomic potential acquires the form of a double well with a 
sizable tunneling rate. The splitting must be performed 
slowly compared to the tunneling rate in order to ensure 
that one atom goes to each site. Note also that, in the 
absence of interaction between the atoms, one cannot 
map a doubly occupied state into two singly occupied 
orthogonal states. Thus, fast non-adiabatic merging and 
separation is not possible within this approach. 



B. Gate using spin-dependent transport 

To realize a faster gate, we propose to adapt the spin- 
dependent lattice transport process [11] to our architec- 
ture by utilizing spin-dependent light shifts in such a way 
that only particular spin combinations are merged. By 
using |F = 2,m F = -2) = |f) and \F = l,m F = -1) = 
\l) and an appropriate laser detuning between the two 
excited state fine structure levels (420.86 nm) the light 
shift for the state \\.) cancels when using o~ polarized 
light. 1 We propose to perform a quantum phase gate in 
five steps (Fig.El). 

We consider an atom in a superposition state of |f) and 
\l) . First, the |f) part of the wave function is moved to an 
empty row. Next, the |f) part of a neighboring qubit is 
moved to overlap with the remaining \\.) of the first atom. 
In this way, a controlled collision between these two spin 
states with an energy Ufi takes place, while the other 
combinations of spin states undergo no interaction. The 
spin state thus acquires a phase shift </>, which attains 
7r after a certain amount of time. Finally, the steps are 
reversed to restore the original configuration. Since each 
spin state resides in its own trapping potential during the 
collision, the two atoms can be rapidly separated after the 
collisional interaction. 

Unlike in the transport discussed in previous sections, 
the phase shift during merging and interaction is acquired 
by only one of the two-qubit spin state combinations. 
Since the phase depends deterministically on the inten- 
sity profile during transport and on the transport time, it 
can be calculated and corrected for. Choosing intensity 



1 Also, using <r+ polarized light, one can cancel the light shift of 
either of the states, but the ratio of light shift and scattering rate 
is worse in this case. 




Figure 4: Phase gate using spin-dependent transport, a) The 
two atoms (blue and red) are in an arbitrary superposition 
between \\) and |t) and are located at neighboring sites after 
the red atom has been transported. Here the tweezer is still 
at the full depth as used for the transport. The lattice site 
to the left of the pair is empty (Fig.fTb. b) While the tweezer 
on the right atom is switched to circular polarization acting 
only on the |t) component, a second tweezer with circular 
polarization is ramped up at the position of the left atom, 
c) The |t) components of both atoms are moved one lattice 
site to the left, d) The polarization of the right tweezer is 
switched to linear, effectively ramping up the potential for the 
|4_) component. The atoms in the central lattice site acquire 
a collisional phase <f>- Subsequently the steps are reversed to 
restore the original configuration. 



step 


amount 


time 


ramp-up/down 


6 


11 /is 


transport 


2(n + l) 


25 /is 


phase gate 


1 


83 /xs 


overall 


- 


(199 + n50) us 



Table I: Time budget for the spin-dependent transport gate 
involving a transport over n lattice sites. 



obtain reasonable coherence times. The dephasing time 
T c = h/AE of the two states is given by their relative 
energy shift, AE = 3/2fi B B = h • 2.1 kHz/mG. It is state 
of the art of active magnetic field stabilization to reach 
a short term stability at sub-milligauss level [3H EZ] . A 
magnetic field noise of 50 /iG yields a coherence time of 
T c = 10 ms and thus leads to a phase error of 10~ 2 within 
the gate time of T g ~ 100 /is. An improvement to 5/xG 
is extremely challenging but could be realizable with ex- 
isting technology [35] and it would lead to an error of 
10~ 3 . 



C. Spin-exchange gate 



and transport profiles where the light shift is symmetri- 
cal with respect to the center of the pulse sequence, a 
robust compensation can also be achieved by applying a 
7r-pulse after half of the sequence to reverse the phase 
evolution and another 7r-pulse at the end to restore the 
correct populations. 

To optimize the speed of the collisional phase gate, 
we switch the optical tweezer to linear polarization after 
the merging to ramp up the tweezer for the other spin 
state and enhance the interaction strength. For a tweezer 
depth of 500 E T the interaction energy is U = h ■ 6 kHz 
and a phase shift of it is acquired in t„ « 83 /is. In Ta- 
ble [I] we summarize the times required for the individual 
steps, obtaining i gate « (200 + n ■ 50)/xs, where n is the 
transport distance in units of lattice sites. In a possible 
implementation, the transport of several atoms could be 
performed in parallel to keep the typical values of n per 
gate low. Although this does not improve the individual 
gate error, it does decrease the effective time per gate and 
thus diminishes the effect of background decoherence ef- 
fects. 

The figure of merit for the speed of a quantum gate 
is the ratio of the gate time to the coherence time 
of the qubit. The hyperfine states which allow spin- 
dependent potentials, such as \F = 1, m F = — 1) and 
\F = 2, m p = — 2) are sensitive to magnetic fields and it 
is therefore crucial to reduce the magnetic field noise to 



When using the clock states \F — l,mp = —1) and 
\F = 2, nip = +1) of 87 Rb as qubit states, coherence 
times of many seconds have been observed working at 
the 'magic' magnetic field of 3.229 G 05]. In the following 
we present a gate proposal based on atoms in these clock 
states. Spin-dependent potentials do not exist for the 
clock states, and we thus need to merge the atoms to the 
same potential. In contrast to the simple merging gate 
of Section IV. A, we propose to transfer the two qubits 
into different vibrational states of the combined poten- 
tial instead of merging them into the ground state. This 
separation in orthogonal spatial wave functions will allow 
more robust merging and separation of the two qubits. 

It can be done by mapping the left well qubit into the 
excited vibrational level of the combined well and the 
right well qubit into the ground vibrational level, while 
preserving the qubit state [12] . This mapping is described 
by: 
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(5) 



where |-) L and |-) R designate the wave function of the left 
and right potential well and |-) , |-) are the ground and 
first excited vibrational states in the right potential well. 

In the combined system the new eigenenergy basis is 



(6) 



now formed by the singlet and the triplet states 

l«) = lt) g ll) c -|l) g lt) c 

|*0> = lt) g ll)e+ll> g lt) c 

l*-l) = \i) e \De 

l«+l) = lt) g |t) e • 

For the v / swap gate, we exploit the differential phase 
evolution between the singlet and the triplet states which 
stems from their different symmetry. As the singlet spin 
state \s) is antisymmetric, the spatial wave function of 
bosonic particles must also be antisymmetric, while the 
symmetric triplet state \t) leads to a symmetric wave 
function. In the antisymmetric wave function of \s), the 
two atoms have essentially zero overlap, such that inter- 
actions play no role. In the symmetric wave function 
of state \to), however, the collisional interaction energy 
has a finite value U eg and a corresponding phase is ob- 
tained. The value of U eg is related to the interaction 
energy U gg of two atoms in the ground state via integrals 
of the appropriate spatial distributions, and for harmonic 
oscillator states U eg — 0.35 • U gg [52] . 

The differential phase evolution between \s) and |£o) 
results in spin-exchange oscillations between \\) \l) c and 
ll)Jt) c : 



tf (i = 0) 

*(* = r swap ) 

*(* = T smp /2) 



+ \to) 

+ e iU ': 

-\to) 
+ i\to] 



- It), ll>e 

t/h |<o) 

- \i) g lt) e 
~\t) \l) t 



ll)Jt) f 



After a time T!3 Wap = ■nfi/U, 



'eg 



(7) 
the states are swapped, 
while after t = T swap /2 an entangling ^/swap gate is re- 
alized. 

For the band mapping process, the two atoms for the 
two-qubit gate are transported to neighboring lattice 
sites [see Fig.[5^)-c)]. The transport tweezer on the left 
atom starts at a depth of Vstart = 400 E r and is linearly 
ramped down to zero in 75 /zs while it is simultaneously 
moved to the right lattice site. We find that the perfor- 
mance of the band mapping is improved if we apply a 
second auxiliary tweezer on the right lattice site, which 
is kept at a constant depth of V aux = 200 E r . This con- 
figuration leads to a mapping of the left atom to the first 
band on the right lattice site [Fig.J5^,)-c)]. 

To optimize the fidelity of the band mapping, we re- 
fine the linear ramp shapes by adding harmonics of the 
ramp time to both the intensity and the position ramp 
of the transport tweezer. We solve the time-dependent 
Schrodinger equation and obtain the single-particle map- 
ping fidelities as the overlap of the calculated wave func- 
tion with the target states. We define the fidelity F of 
the merging process as the product of the single-particle 
mapping fidelities of the two atoms. Fig.[HH)-g) shows 
the resulting ramps and wave functions for an optimiza- 
tion including harmonics up to 1 5th order corresponding 



step 


amount 


time 


ramp-up/down 


2 


11 /iS 


transport 


2-ra 


25 /is 


phase gate 


1 


125 /zs 


merge/split 


2 


75 /zs 


overall 


- 


(297 + n ■ 50) /zs 



Table II: Time budget for the exchange gate. 



to a spatial period of 0.13 ai at , 
infidelity of 1 - F = 2 • 10" 4 . 



for which we obtain an 



A similar analysis for double- well lattices resulted in a 
saturation of the merging error of the order of 10 -2 |19j . 
which is significantly higher than our results. We at- 
tribute this to the difference in available control param- 
eters 2 . We have also included interactions into the sim- 
ulations but the transfer is so fast that the effect on the 
overall band mapping fidelity is negligible. It is impor- 
tant to calculate the accumulated interaction phase dur- 
ing the band mapping as it contributes to the total gate 
phase. In the case illustrated in Fig. [3] the accumulated 
interaction phase during the mapping was 0.15 rad. 

After the band mapping, we wait for a time T swap /2 
to perform the ^/swap gate before reversing the band 
mapping process. For comparison with the gate proposal 
discussed in Sec. IV B we introduce an 11 /is ramp up of 
the tweezer to 500 E T) which results in a reduction of the 
swap time to T swap /2 = 125 /xs. 

In summary, the exchange gate requires four ramps 
of the tweezer power (each ~ 11 /xs) , the transport 
(25/x/site), the merging and splitting of the wells 
(150 /xs), and finally the phase gate (125 /xs) (see Tab. III]). 
In total this gives a gate time i gatc ~ (300 + n ■ 50) /zs, 
where n is the distance to the interaction site. The gate 
time is therefore more than 10 4 times shorter than the 
decoherence time of the clock states. 



V. Discussion of errors 

In the following we discuss the contribution of error 
and decoherence sources other than the vibrational exci- 
tations considered so far. 



2 In Ref. |19j the control parameters are the overall power of the 
double well lattice, the ratio of horizontal to vertical power, and 
the phase between the two, which indirectly determine the posi- 
tion and relative depths of the two sides of the double well. 
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Figure 5: Band mapping for the exchange gate. a)-c) 
Schematics of the potential formed by the lattice and the two 
tweezers at three different times. While the transport tweezer 
is moved from the left to the right lattice site and simultane- 
ously ramped down, the atom in the left well is mapped to 
the first excited state in the right site, d) and e) Depth and 
position of the transport tweezer as a function of time for the 
optimized ramp shapes. The auxiliary tweezer is kept at a 
fixed depth and position (not shown), f) and g) Density pro- 
file versus time for the transported atom f) and the stationary 
atom g). At the end of the sequence the atom starting in the 
left well is mapped to the first excited state and has a wave 
function with a node, while the atom starting in the right well 
is unchanged. 



A. Spontaneous scattering 

A fundamental source of error is the spontaneous scat- 
tering due to the off-resonant absorption of lattice and 
optical tweezer photons. Due to the large detuning, the 
absorption of lattice photons can be neglected. With a 
scattering rate of T sc w 4 • 10 -2 Hz for all three lattice 
axes, the probability to scatter a photon during one gate 
time is < 10~ 6 . 

For the optical tweezer, the spontaneous emission can 
be suppressed by choosing a large detuning. A detuning 
of 10 nm requires a laser power of 1 mW to reach the 
500 E T potential depth. This yields a scattering rate of 
T sc = 0.1 Hz and a scattering probability during one gate 
time T g ~ 300 ^s of P sc = r sc T 9 = 3 • 10" 5 . 

In the case of spin-dependent transport, the detun- 
ing is much smaller since the wavelength of the optical 
tweezer has to lie between the two fine structure lines. At 
the chosen waist radius, the potential depth of 500 E t is 
reached at a power of 70 /xW yielding a spontaneous scat- 
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Figure 6: Band-mapping infidelity as a function of intensity 
noise of position errors. For the transport tweezer we intro- 
duce a positioning error and a scaling error in the amplitude. 
The contour lines show the multiples of an error of 10 -3 . 



tering rate of r sc = 1.5 Hz and a scattering probability 
of P sc = 4 • 10~° during one transport time. 



B. Pointing stability and intensity noise 

We now discuss the influence of experimental imper- 
fections on the gate fidelity. For the present proposal the 
limitations are the pointing error and intensity fluctua- 
tions of the tweezer. In order to investigate these effects 
quantitatively, we assume only low-frequency noise due 
to shot-to-shot fluctuations. We calculated the gate in- 
fidelity in the presence of a constant offset in the point- 
ing position and of a reduction or increase of the opti- 
cal tweezer intensity by a constant factor. In particular, 
we included these errors into the simulation of the band 
mapping process discussed in Sec. IV C. As can be seen 
in Fig.|6j both parameters have to be very well controlled 
in order to avoid a significant reduction of the fidelity. In 
order to keep the gate infidelity below 10 -3 , our scheme 
requires a relative intensity stability of 10 -3 and pointing 
precision in the nm range. This intensity stability is defi- 
nitely manageable whereas the required pointing stability 
is beyond the current state-of-the-art experiments, which 
achieve a pointing error of 50 nm [T7j, where full active 
pointing stabilization was, however, not implemented. It 
might be possible to find conditions under which the per- 
formance is less sensitive to errors in the input controls 
as it was done for transport of harmonic traps [50] . 

Intensity noise causes additional dephasing in the spin- 
dependent two-qubit gate of Section V.B, since the depth 
of 500 E T = h ■ 1 MHz implies a phase accumulation of 
several tens of 27rrad during the interaction time. To 
keep the absolute dephasing error below 2ir ■ 10~ 3 rad, 
the relative intensity noise has to be below 10~ 5 . This 
is technically challenging but is possible as demonstrated 



in the LIGO project [51], where the intensity noise is 
reduced to a fractional level of 10 -8 over the relevant 
spectral band. 



C. Error threshold for fault-tolerance 

We have argued that combined transport and gate er- 
rors at the 10 -3 level can be reached with our proposal. 
With errors of this magnitude, error correcting codes may 
be applied to reach scalable fault tolerance, but only with 
a vast overhead in number of physical qubits and oper- 
ations [52] . An alternative approach to fault tolerance 
uses stabilizer codes in a scheme where errors are not 
corrected, but error syndromes are forwarded to higher 
algorithmic levels, and here gate errors at the few per- 
cent level can be tolerated |53) . Fault tolerant quantum 
computing is a very active field of research, and very re- 
cent simulations show tolerance to percent level errors 
for a two-dimensional qubit structure on a lattice with 
only nearest neighbor interactions |54j . With our ability 
to move atoms around with low error, a simpler scheme 
with the same error threshold may well be possible. 



VI. Conclusion 

In conclusion, we have proposed a complete architec- 
ture for scalable quantum computing based on transport 
of atoms in an optical lattice by movable optical tweez- 
ers. We have presented explicit calculations of the errors 
occurring in the different processes of our proposals, and 
we have numerically optimized non-adiabatic ramps of 
the tweezer potentials to a vibrational excitation error 
below 10 -3 . Our analyses show that ultracold atoms 
in optical lattices are indeed a promising candidate for 
fault-tolerant scalable quantum computing. 
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